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ABSTRACT 

We present limits on planetary companions to pulsating white dwarf stars. A subset 
of these stars exhibit extreme stability in the period and phase of some of their pulsation 
modes; a planet can be detected around such a star by searching for periodic variations 
in the arrival time of these pulsations. We present limits on companions greater than 
a few Jupiter masses around a sample of 15 white dwarf stars as part of an on-going 
survey. One star shows a variation in arrival time consistent with a 2 Mj planet in a 4.5 
year orbit. We discuss other possible explanations for the observed signal and conclude 
that a planet is the most plausible explanation based on the data available. 



Subject headings: planetary systems — wliite dwarfs 



1. Introduction 

All main-sequence stars with mass less than 
about 8 Mq will end their lives as white dwarf 
stars (WDs). As such, WDs are a fossil record 
of star formation in the Galaxy from earliest 
times to just a few hundred million years ago. 
WDs offer a window into the ultimate fate 
of planetary systems, including our own solar 
system, and whether planets can survive the 
final stages of stellar evolution. The proper- 



ties of a WD are relatively insensitive to the 
mass of the progenitor: the mass distribu- 
tion of isolated WDs is narrowly distributed 
in a peak 0.1 Mp) wid e around a mean mass of 



0.59 M0 (jKepler et al.ll2007l l . Surveys of WDs 
can therefore search a wide range of main se- 
quence progenitor masses and the low lumi- 
nosity of a WD means any companion planets 
can be potentially followed up with current di- 
rect detection technology in the mid-infrared. 
Livio &: Soker ( 19841 ) considered the fate of 
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a planet engulfed in the envelope of a red gi- 
ant star. They determined that below a cer- 
tain mass the planet would be evaporated and 
destroyed, but larger objects would accrete 
material and spiral in toward the stellar core. 
They predicted that the end state of these sys- 
tems would be a tight binary consisting of a 
white dwarf and a brown dwarf and suggested 
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this mechanism might explain the origin of 
cataclysmic variable (CV) systems. 

Planets that are not engulfed, and suffi- 
ciently far from the stellar surface that tidal 
drag is small, will drift outwards to co nserve 
angu lar momentum (as des c ribed in [Jeans 



1924 ). iDuncan &: Lissaueii (|l998l ) investi- 



gated the stability of the outer solar sys- 
tem when the sun undergoes mass loss as 
a red giant. They found that the system 
was stable on timescales of at least 10 bil- 
lion years for reasonable amounts of mass 
loss, but for larger amounts typical of WDs 
formed from more massive stars, the planets' 
orbits be came unstable on tini e scales of < 10^ 
years. iDebes fc Sigurdssonl (|2002l ) looked 



what would happen if the orbits of two plan- 
ets became unstable and determined that if 
orbits crossed, the most likely result was that 
one planet would be scattered into a shorter 
period orbit, while the other would be boosted 
into a longer orbit or ejected from the system. 
For planets scattered inward, the extreme flux 
from a newborn white dwar f would strip the 
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outer atmospheric layers. 

hooi ) estimated that a 2Mj planet 1.8 AU 
from a young white dwarf would lose half its 
mass in this manner. The picture drawn by 
this brief survey of theory is a population of 
planets in long period orbits around WDs, 
with a number of objects scattered closer to 
the star, or in very tight binaries. 

Searches for sub-stellar companions to 
WDs have concentrated on exploiting the 
lower contrast between star and c ompan- 
ion, especially in the infrared (e.g 



1983: Burleigh et al.1 2002: Farihi et al 



'robst 



2005 



report the detection of a planetary mass com- 
panion in a 1.7 AU orbit around an extreme 
horizontal branch sdB star. 

Pulsating WDs allow the possibility of 
searching for the presence of companion plan- 
ets as changes in the observed arrival time of 
the pulsations. Hydrogen atmosphere (DA) 
WDs pulsate in an instability strip approx- 
imately 1200 K wide near 12,000 K and are 
known as DAVs (or ZZ Ceti stars). The pul- 
sations are non-radial g-modes with periods 
of order 100-1500s and amplitudes of a few 
percent. 

The pulsation properties of DAVs vary with 
temperature. Those near the hot end of the 
strip tend to have a smaller number of shorter 
period, lower amplitude modes with sinu- 
soidal lightcurves and are known as liDAVs, 
while those nearer the red edge show more, 
larger amplitude modes and asymmetric pulse 
shapes. The change in pulsation properties is 
most likely due to increasing depth of the 
convection laye r near the surface of the star 
(|Brickhil]||l9"8a . and subsequent articles). 



Stover et al. I (|l980l l first noted that modes 



on hDAVs often exhibited an impressive sta- 
bility in the period and phase of pulsation. 
Kepler et al. ( 20051 ) measured the rate of pe- 
riod change, P, of one mode in the hDAV 
G117-B15A to be 3.5 7(82) x 10~^^, while 
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mode in R548 



(|2003l l constrained 
to < 5 5(1.9 ) X 10" 



Friedrich et al.ll200fil : lDebes et al.ll200fil : lMullallv ^ 
20071). Although a couple of brown dwarf 
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stars have been found with this approach 
(IZuckerman &: Becklinlll992l : lFarihi &: Christopher 
20041 ). as yet no di rect detections of pl anets 
have been claimed. Silvotti et al. ( 20071 ). us- 
ing the same timing method discussed here. 



P of a 
"^^ and 

O'Donoghue &: Warneil ( 1987 ) constrained P 
of L19-2 to < 3.0 X 10"^^. Together with 
pulsars, these objects are the most stable as- 
trophysical clocks known. 

Measurements of P require datasets of be- 
.ween 10 and 30 years, but the investment in 
time yields a suitable scientific reward. Mea- 
surement of P provides a rare opportunity 
to directly test models of struct ure and com- 



posit ion of the core of a star (iKepler et al 



20051 ). constrain the current rate of change of 



the gravitational constant (jBenvenuto et al 
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2004 ) . as well as provide useful constraints on 
the mass of the hypothesiz ed axion or other 
super-symmetric particles (llsern et al. I1992I : 
Corsico et allboOll : iBischoff-Kim et al.l[2007h 
If a planet is in orbit around a star, the 
star's distance from the Sun will change pe- 
riodically as it orbits the center of mass of 
the planetary system. If the star is a stable 
pulsator like a hDAV, this will cause a peri- 
odic change in the observed arrival time of the 
otherwise stable pulsations compared to that 
expected based on the assumption of a con- 
stant period. The change in arrival time, r, is 
given by 

apmpsini 

where Up is the semi-major orbital axis of the 
planet, rup is the planet mass, M^, is the mass 
of the white dwarf, c is the speed of light and 
i is the inclination of the orbit to the line of 
sight. In common with astrometric methods, 
the sensitivity increases with the orbital sep- 
aration, making long period planets easier to 
detect given data sets with sufficiently long 
baselines. 

In 2003 we commenced a pilot survey of 
a small number of DAVs in the hope of de- 
tecting the signal of a companion planet. We 
present here a progress report of the first 3-4 
years of observations on 12 objects, as well as 
presenting limits around 3 more objects based 
partly on archival data stretching as far back 
as 1970. For one object we find a signal con- 
sistent with a planetary companion. Further 
observations are necessary to confirm the na- 
ture of this system. For our other objects, we 
can constrain the presence of planets down to 
a few Jupiter masses at 5 AU, with more strin- 
gent limits for stars with archival data. 

2. Our Survey 



didate DAVs with spectra taken by the Sloan 



Digit al Sky Survey (jAdelman-McCarthv et al 



20061') . A fol low up survey b y 



Mukadam et al 



tooi ) and iMullallv et aD (|2005l ^ confirmed 
46 of these candidates to be pulsators. We 
selected our targets from these two papers 
as well as earlier kn own DAVs published 
in the literature (s ee iBergeron et al.l Il995l : 
Fontaine et alll2003l l. 



.Kleinman et al.l (|2004l ) and lEisenstein et al 
(j2006l ) published a large number (~ 600) 



can- 



The ideal DAV for this survey would ex- 
hibit a number of isolated, relatively low am- 
plitude (0.5-2%) modes. Multiplet, or oth- 
erwise closely-spaced (< 1 s) modes are diffi- 
cult to resolve in single-site data and inter- 
ference between the unresolved modes makes 
accurately measuring the phase difficult. We 
selected a sample of 15 stars brighter than 
]^gth jj^g^g fQj. long-term study. With the ex- 
ception of R548, which has a well studied 
double mode, we chose stars with one iso- 
lated mode with amplitude > 0.5%. Only 
one star, SDSS J221458.37-002511.7 has two 
modes suitable for study. 

We have monitored this sample of stars for 
4 years us ing the Argos Prime Focus CCD pho- 
tometer (jNather fc MukadamI |2004| ) on the 
2.1m Otto Struve Telescope at McDonald Ob- 
servatory. We list the observed objects in Ta- 
ble dl along with the period and P of the 
analyzed modes. We observed each object 
with exposure times of 5-15s for periods of 4-8 
hours per night. The exposure times are cho- 
sen to be very much shorter than the Nyquist 
frequency of the shortest period mode ob- 
served on the star (> 100s), and the observing 
time to sample many consecutive cycles of the 
pulsation. 

We re duce our data in th e manner de- 
scribed in iMullallv et (|200,5l ) with one im- 
provement. Argos suffers from a fiuctuating 
bias level but does not have an overscan re- 
gion. To account for this we measure the bias 
from a dead column and subtract this value 
from both our science and dark frames be- 
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fore flatfielding. This is clearly not ideal, but 
the best approach to measuring the bias avail- 
able, given that the level can vary by up to 5 
DN/pixel on timescales shorter than the expo- 
sure time. 

We perform weighted aperture photometry 
with a variety of apertures using the IRAF 
package wphot, choosing the aperture that 
gives the best signal-to-noise by eye. We di- 
vide the light curve by a combination of one 
or more reference stars, remove points affected 
by cloud, fit a second order polynomial to re- 
move the long term trend caused by differen- 
tial extinction, and correct our timings for the 
motion of the Earth around the bary center of 
the so lar system using the method of lStumpfl 
( IQSd ). accounting for all UTC leapseconds up 
to and including January 2006. 

We combine all the data on a star in a given 
month for analysis, typically 8-16 hours. We 
first compare the alias pattern of the peaks 
in the Fourier transform (FT) with a window 
function to identify closely spaced modes and 
multiplets. It is more difficult to measure the 
phase of closely spaced modes (< 70 //Hz) be- 
cause interference between the unresolved pe- 
riodicities requires significantly more data to 
resolve, so we focus only on isolated modes. 

Having selected a mode for analysis, we 
first measure the period by fitting a sine curve 
using the Marquant- Levenberg non- linear 
least squares technique ( Bevington 19691 ). We 
attempt to obtain between two and four ac- 
curate timings on each star per year. As we 
accumulate data, we re-measure the period 
using the entire dataset, before measuring the 
phase of that period in each month's data 
using a least squares fit. We compare the ob- 
served phase to that expected based on the 
assumption of a constant period and plot the 
result in an 0-C diagram. We also check that 
the amplitude of a mode is stable from month 
to month. Varying amplitudes are a symp- 
tom of either unresolved companion periods. 



or an instability in the pulsation mechanism. 
None of the modes discussed here displayed 
any amplitude variability inconsistent with 
observational error. 

3. Results 

3.1. GD66 

GD66 (WD0517-F307) is an 11,980 K, 
log q =8.05, 0.64 Mq hDA 



2004) 



(Bereeron et al. 


of 15.6 ( 


Eescen 



19681 ) corresponding to a dista nce of about 
51pc (jMullallv fc de Grafj [2OO5I 1 . The FT is 
dominated by a single mode at 302 s, triplets 
of modes separated by ~6.4 /iHz at 271 and 
198 seconds and a handful of other lower am- 
plitude modes. There are also some combina- 
tion and harmonic peaks present. A sample 
FT is shown in Figured) 

The presence of closely spaced peaks at 271 
and 198 s makes it more difficult to accurately 
measure their phase and our analysis concen- 
trates on the 302 s mode. We show an 0-C di- 
agram for the arrival times of the 302 s mode 
in Figure [2j The curvature in this diagram 
is unmistakable. Instabilities in the pulsation 
modes of DAVs often manifest as variations 
in the amplitude of pulsation. In Figure [3] 
we plot the amplitude of the 302s mode as a 
function of time and find the amplitude varies 
between 1.1 and 1.2%. We can reproduce vari- 
ations of this magnitude by small changes in 
our reduction method, and conclude that the 
amplitude is stable within our ability to mea- 
sure it. By fitting a sine curve to the 0-C 
diagram we find a period of 4.52 years and an 
amplitude of 3.84 s. 

The data collected to date is also con- 
sistent with a parabola. 0-C diagrams of 
hDAVs are expected to show parabolic be- 
havior as the cooling of the star produces a 
monoton ic increase in the period of pulsa- 

However, based 



tion (see lKepler et al.lll991 

on observations of other DAVs and mod- 
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els of white dwarf in teriors ( Bradley 19981 : 
Benvenuto et al. 20041 ) we expect the cooling 
to cause a -Pcooi ~ 10~^^. If we fit a parabola 
to our data we find a P = 1.347(95) x 10" 
three orders of magnitude larger than ex- 
pected from cooling alone. 

The tangential motion of the star with 
respect to the line of sight also causes a 
parabolic curvature in the 0-C diagram. As 
the star moves linearly in space perpendicular 
to our line of sight, its distance to us ch anges 
parabolically ( Shklovskii 197d: Paidosz 1995i) 



dWinget et al.l Il990l l. However, analysis of 
other modes on th e same star failed to r epro- 
duce this behavior ( Kleinman et al. 199^1 and 



The USNO-Bl catalogue (jMonet et al 



20031 ) 



quotes a proper motion of 131.6(5.0) mas/yr 
corresponding to a Pp^ of 6.4 x 10~^^, again 
too small to explain the observed data. 

Apparently periodic signals in 0-C dia- 
grams can be caused by random jitter or drift 
in the period of the pulsator. A likelihood 
statistic, C, that a given data set was caused 
by different combinations of observational er- 
ror, period ji tter and drift can be calculated 
according to Koen ( 20061 ). We expand his 
methodology to determine the likelihood that 
the data shown in Figure[2]is the signature of a 
companion, or the result of stochastic changes 
in the pulsation period. We first calculate C 
for a model of the data that seeks to explain 
the data by invoking jitter or drift in the pe- 
riod and find values for log£ of -15.1673 and 
-15.0918 respectively. A model including both 
jitter and drift gives a similar value. Next, 
we calculate the likelihood that the residuals 
of the sine fit can be explained by observa- 
tional error alone, and find a value of log £ of 
-11.5549. This strongly disfavors the hypoth- 
esis that the observations can be explained by 
small random changes in the pulsation period. 

Time series observations of another DAV, 
G29-38, showed a variation in the phase of 
one mode over 3 months consistent with an 
O.5M0 object in an eccentric 109 day orbit, 
but a change in amplitude the following year 
made the mode unreliable as an accurate clock 



near-infrared imaging ( Kuchner et al. 19981 : 



Debes et al.l |2005| ) did not detect any sub- 
stellar companions. It is possible that the 
same internal effect that mimicked a compan- 
ion to G29-38 is also present in GD66 albeit 
with a much smaller amplitude and consider- 
ably longer period. 

If we assume the curvature is caused by 
a planet in a circular orbit, the best fit pe- 
riod is 4.52(21) years and we can use Ke- 
pler's laws to determine an orbital separa- 
tion, Up = 2.356(81) AU. The amplitude of 
the sine curve, r = 3.84(32) s is related to 
the semi-major axis of the star's orbit, a^:, by 
Tc = a=K sin i where i is the inclination of the 
orbit to the line of sight and c is the speed 
of light. The mass of the planet, nip, is equal 
to {M^a:^) / ap, where M* is the mass of the 
star. Using these two equations we find an 
mpsini of 2.11(14) Mj. From our best fit cir- 
cular orbit, we predict we will obtain observa- 
tions spanning an entire orbit in early 2008. 

3.2. Other stars 

GD66 is the only star in our sample that 
shows strong evidence for a planetary com- 
panion. However, we can place interesting 
limits on the presence of planets around the 
other stars. 

The fundamental limit on our ability to de- 
tect planets is set by the scatter in the 0-C 
diagram. Other factors which affect this limit 
include the timespan and sampling pattern of 
the data. We perform a Monte Carlo analy- 
sis to estimate the region of the mass-orbital 
separation plane in which planets are actually 
detectable around each star based on the data 
available. We randomly choose a planet mass, 
orbital separation, eccentricity, and other or- 
bital parameters and calculate the effect this 
planet would induce on the 0-C diagram of a 
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0.59 Mq white dwarf. We then sampled this 
0-C diagram with the same observing pattern 
and error bars as our actual data for each star 
and fit the resulting 0-C diagram with a sine 
curve and a parabola. If either the ampli- 
tude of the sine curve or the curvature of the 
parabola were measured with 2>a confidence, 
the hypothetical planet was determined to be 
detected. We repeated this process 10^ times 
for each star and drew a shaded relief map in- 
dicating the percentage of the time a planet 
with a given mass and orbital separation was 
detected with either technique to Scr, with 
dark shades indicating near 100% detection 
efficiency, and white indicating regions where 
planets were unlikely to be detected. We show 
the 0-C diagram and the relief map for each 
star in Figures 3-16. As can be seen in the fig- 
ures, the annual sampling pattern means that 
planets with orbital periods of 1 (Earth) year 
are difficult to detect, resulting in the weak 
limits for planets at separations of slightly less 
than 1 AU. 

3.3. Notes on Individual Stars 

Gir/-B15A.— Also known as WD 0921-H354. 
30 yea rs of archival data comes from lKepler et aP 
(|2005l l. although some of the more recent 
data in that work was taken in conjunction 
with this project. Where data did not come 
from our observations we used the 0-C value 
quoted in their Table 1. With this timebase, 
the curvature caused by the change in period 
due to cooling becomes evident. This cooling 
effect is removed from the data before per- 
forming the Monte Carlo simulation. The lim- 
its on long period planets placed around this 
star are among the best constraints placed 
around any stellar object by any technique. 
See Figure m 

G185-32.~ Also known as WD1935+279. 
The point from the early 90 's comes from 
archival data from the Whole Ea rth Telescope 
(XcovS. ICastanheira et al.l 12004 ). This extra 



point gives a long baseline, but the poor cov- 
erage in our survey reduces the sensitivity. 
See Figure m 

R548.— Also known as ZZ Ceti and WD0133- 

116. T he entire data set comes from lMukadam et al, 

(|2003l ). See Figure [3 

SDSS JOlllOO. 63+001807.2.— At g =18.6*'^ 
magnitude, this is our faintest target and cor- 
respondingly has our weakest limits. Also, 
because it could only be observed under the 
best conditions the data coverage is quite low. 
The apparently impressive curvature in this 
0-C diagram in entirely due to the last data 
point and should be treated with considerable 
skepticism. See Figure [TOl 

SDSS JI35459. 89+010819. 3.— This bright 
(16.4*^ magnitude) star has a baseline stretch- 
ing back to early 2003 and some of the highest 
accuracy time measurements, and as a result 
has the best limits on planets for stars with- 
out archival data. We could have detected 
a Jupiter mass planet at 5AU had one been 
present. See Figure [T4l 

SDSS J221458.37-002511.7.— This is the 
only star in the sample for which reasonable 
0-C diagrams were obtained for two modes. 
The 0-C diagram shown (see Figure [T7|) is the 
weighted sum of the 0-C values for these two 
individual modes. 

4. Discussion 

In this pilot study we find a signal con- 
sistent with a companion planet around one 
star in a sample of only 15. Although further 
observations will be necessary to conclusively 
identify the cause, it augers well for the fu- 
ture potential of this and other white dwarf 
planet searches. According to the theoretical 
arguments discussed in 311 the distribution of 
planets around WDs can be expected to be 
weighted toward planets in long period orbits. 
The population of hot Jupiters around the 
main sequence progenitors will likely be de- 
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stroyed, while more distant planets will drift 
outward with stellar mass loss. Evidence of 
a planet in a relatively short orbit encour- 
ages us to continue monitoring for planets 
with greater orbital separations. For the other 
stars in our sample, we can rule out the pres- 
ence of planets down to a few Jupiter masses 
at 5AU. With more data, we will extend our 
limits beyond 10 AU into the regime where we 
expect planets to be most frequent. 

Our search spans a broad range of progeni- 
tor masses. To estimate the progenitor masses 
we compare the spe ctroscopically in e asured 
Teff and log g froni Eiseristein et al. ( 20061 ) 
and iBergerori et al.l (120041') to W E > inte rior 
models from Holberg &: Bergeron ( 20061 ) to 
obtain a WD mass. We then used the initial fi- 



nalmass rel ations (IFMRs) of I Williams et al 



(|2004l ) and iFerrario et aD (j2005|) to calcu- 
late the progenitor mass, taking the weighted 
mean of the two relations as the best value, 
and adding the difference between the two 
methods in quadrature to the uncertainty. 
We present the masses in Table [2l Although 
there is still considerable uncertainty in this 
relation, the results give some indication of 
the type of star that created the WD. The 
progenitor masses correspond to a range of 
spectral type frorn appro ximately B6 to F9 



(jHabets Sz Heintzd Il98ll ). a range that is 



largely complementary to the radial velocity 
method. 

5. Conclusion 

We present our results on an on-going sur- 
vey for planets around 15 pulsating white 
dwarf stars. Our survey data spans 3-4 years 
with archival data on some stars stretching 
back to 1970. We are already sensitive to 
planets down to a few Jupiter masses at dis- 
tances of approximately 5AU. For one star, 
we observe a curvature in the 0-C diagram 
consistent with a planet in a 4.5 year orbit. 
Further observations are necessary to span a 



full orbit of this candidate object. If con- 
firmed, this will be the first planet discovered 
around a WD, and together with the planet 
discovered around an sdB, suggests that plan- 
ets regularly survive the death of their parent 
star, and that WDs will be fruitful targets for 
planet searches. 

This work was supported by a grant from 
the NASA Origins program, NAG5-13094 and 
is performed in part under contract with the 
Jet Propulsion Laboratory (JPL) funded by 
NASA through the Michelson Fellowship Pro- 
gram. JPL is managed for NASA by the Cal- 
ifornia Institute of Technology. 
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Table 1 

Modes used to construct 0-C diacrams 





Star 


Period 

(sec) 


Amplitude 

% 


To 
(bjd) 


P 








G117- 


-B15A. . . 


215.1973888(12) 


1.9 


2442397.9194943(28) 


-1.07(49) 


X 


10~ 


-13 


G185- 


-32 


370.2203552(55) 


0.1 


2453589.6557652(39) 


-0.5(1.0) 


X 


10" 


-13 


G238- 


-53 


122.1733598(38) 


0.2 


2453168.6334567(35) 


-5.7(2.4) 


X 


10" 


-13 


GD244 


202.9735113(40) 


0.4 


2452884.8712580(31) 


0.2(2.8) 


X 


10" 


-13 


GD66 




302.7652959(21) 


1.2 


2452938.8846146(28) 


1.347(95) 


X 


10" 


-12 


R548 




212.76842927(51) 


0.4 


2446679.833986 


1.2(4.0) 


X 


10" 


-15 


SDSS 


J001836.11+003151.1 


257.777859(13) 


0.6 


2452962.6358455(41) 


9.4(9.2) 


X 


10" 


-13 


SDSS 


JOlllOO.63+001807.2 


292.9445269(90) 


1.9 


2452963.7174455(44) 


3.87(43) 


X 


10" 


-12 


SDSS 


J021406. 78-082318.4 


262.277793(11) 


0.6 


2452941.7929412(37) 


-1.5(7.5) 


X 


10" 


-13 


SDSS 


J091312.74-H403628.8 


172.605159(15) 


0.3 


2453024.8275265(47) 


9.6(9.8) 


X 


10" 


-13 


SDSS 


J101548.01-H030648.4 


254.9184503(56) 


0.7 


2453065.6152116(41) 


7.2(3.6) 


X 


10" 


-13 


SDSS 


J135459.88-h010819.3 


198.3077098(14) 


0.6 


2452665.9507137(33) 


-5.3(7.8) 


X 


10" 


-14 


SDSS 


J135531.03-h545404.5 


323.9518703(69) 


2.2 


2453082.8582407(39) 


1.39(47) 


X 


10" 


-12 


SDSS 


J172428.42+583539.0 


335.536871(14) 


0.6 


2453139.8477241(37) 


1.23(85) 


X 


10" 


-12 


SDSS 


J221458. 37-002511. 7 


195.1406388(64) 


0.4 


2452821.8513218(35) 


6.2(3.6) 


X 


10" 


-13 


SDSS 


J221458. 37-002511. 7 


255.1524057(30) 


1.3 


2452821.8521749(35) 


1.7(2.1) 


X 


10" 


-13 



Note. — To is the time of the arbit arily defined zero t h pul se and is given in units of barycentric corrected 
juhan day. Data on R548 comes from ^Mukada m et al.l (|2003l ) who do not provide a value for uncertainty in 
Tq. Except for GD66, we do not claim statistical significance for the measurement of P for any star. 
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Table 2 
Stellar Parameters 



Star 


Magnitude 


Tcff 
(K) 


log 5 


Initial Mass 
(Mo) 


Final Mass 
(Mo) 


G117-B15A 


15.7 


11630 


7.98 


1.69(51) 


0.595(29) 


G185-32 


13.0 


12130 


8.05 


2.10(44) 


0.638(32) 


G238-53 


15.5 


11885 


7.91 


1.36(55) 


0.562(25) 


GD244 


15.6 


11645 


8.01 


1.85(49) 


0.611(32) 


GD66 


15.6 


11980 


8.05 


2.10(44) 


0.638(32) 


R548 


14.2 


11894 


7.97 


1.65(51) 


0.591(28) 


SDSS J001836.11+003151.1 


17.4 


11696 


7.93 


1.45(54) 


0.571(25) 


SDSS JOlllOO.63+001807.2 


18.8 


11507 


8.26 


3.33(37) 


0.769(37) 


SDSS J021406.78-082318.4 


17.9 


11565 


7.92 


1.40(55) 


0.566(26) 


SDSS J091312.74+403628.8 


17.6 


11677 


7.87 


1.17(60) 


0.542(25) 


SDSS J101548.01+030648.4 


15.7 


11584 


8.14 


2.62(37) 


0.693(32) 


SDSS J135459.89+010819.3 


16.4 


11658 


8.01 


1.85(49) 


0.611(32) 


SDSS J135531.03+545404.5 


18.6 


11576 


7.95 


1.61(61) 


0.580(49) 


SDSS J172428.42+583539.0 


17.5 


11544 


7.90 


1.31(57) 


0.556(25) 


SDSS .]22;lir)8. 37-002511. 7 


17.9 


11439 


8.33 


3.75(37) 


0.811(32) 



Note. — Sloan magnitudes are in the g filter, V magnitudes for the other stars are 
taken from Simbad. The sources of the temperatures, gravities and masses are discussed 
in the text. 
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Fig. 1. — Sample Fourier transform of GD66 from a single 6 hour run. The larger amplitude modes 
are labeled with their periods. The peaks at 271 and 198 seconds are composed of triplets of closely 
spaced modes separated by approximately 6.4 fiHz that are not resolved in this data 




Fig. 2. — 0-C diagram of the 302s mode of GD66. The solid line is a sinusoidal fit to the data. 
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f3.ps 

Fig. 3. — Amplitude of the 302s mode of GD66 as a function of time. The error bars are the formal 
errors of a non-linear least squares fit, the systematic error is approximately 0.1%. An unstable 
amplitude would indicate that the observed phase variations are due to some process internal to 
the star, however the amplitude is stable within our ability to measure it 




1980 1990 2000 0.01 0.1 1 10 100 

Time (years) Orbited Separation (AU) 

f4.ps 

Fig. 4. — Left: 0-C diagram of G117-B15A. Each point represents the phase of pulsation based on 
a single night of data. The solid line is a parabolic fit to the data. Right: Relief map of the region 
of parameter space around G117-B15A where we are sensitive to planets. Dark regions indicate a 
high probability a planet would have been detected. The large white circles indicate the location 
of Jupiter and Saturn, while the small circles mark the postions of known extra-solar planets. 
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Fig. 5. — Same as Figure HI except for a different star, G185-32. In this 0-C diagram, the points 
indicate the phase of pulsation measured from the combination of data taken on this star over a 
month or more as described in the text. 




Fig. 6.— G238-53 
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Fig. 7.— GD244 




Fig. 8.— R548 



14 





2004 S004.5 2005 2005.5 2006 2006.5 2007 
Time (years) 

f9 



0.01 0.1 1 10 

Orbital Separation (AU) 



100 



Fig. 9.— SDSS J001836.11+003151.1 




Fig. 10.— SDSS JOlllOO.63+001807.2 
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Fig. 11.— SDSS J021406.78-082318.4 




Fig. 12.— SDSS J091312.74+403628.8 
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Fig. 13.— SDSS J101548.01+030648.4 




Fig. 14.— SDSS J135459.88+010819.3 
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Fig. 15.— SDSS J135531.03+545404.5 




Fig. 16.— SDSS J172428.42+583539.0 
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